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A theoretical investigation is made of acoustic wave propagation in a periodically stubbed waveg- 
uide. In general the waveguide segments and stubs are made of different materials. The acoustic 
wave in such a system has two independent polarizations: out-of-plane and in-plane modes. The 
band structure and transmission spectrum is studied for diverse geometries using a simple and ef- 
ficient version of the transfer-matrix method. For the same material between the waveguide and 
symmetric stubs the width of some gaps can change, upon varying the stub length or width, by 
more than one order of magnitude. A further modulation can be achieved for different material 
between the stubs and the main waveguide or if the stubs are asymmetric. The gaps in the band 
structure of an infinitely long system correspond to those in the transmission spectrum of the same 
system but with finite number n of units. For n finite i) there exist pseudogaps that gradually turn 
into complete gaps with increasing n, and ii) the introduction of defects gives rise to states in the 
gaps and leads to transmission resonances. 
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I. INTRODUCTION 



The term "band-gap engineering" is well known from decades of research in semiconductors. The recently discovered 
periodic dielectric structures, which exhibit a photonic band gap (PBG) analogous to the electronic band gap in 
semiconductors, have attracted considerable attention due to many interesting phenomena and potential applications 
emerging from them, such as the control of spontaneous emission of radiation, zero-threshold lasing, and the sharp 
bending of light jjj . 

It did not take long before the study of PBG materials, involving light waves, led to analogous studies in other 
systems involving elastic/acoustic waves, e.g., the phononic crystals j| or other periodic acoustic composites ||. The 
phononic crystals ^ have drawn comparatively greater attention, both theoretically [Q and experimentally H[|. 
In analogy with PBG crystals, the emphasis in phononic crystals has been on the occurrence of complete acoustic 
gaps within which the sound, vibrations, and phonons are all forbidden. This is of interest for applications such as 
ultrasonic filters, noise control, and improvement in the design of transducers, as well as for fundamental physics 
concerned with the Anderson localization of sound and vibrations . 

The purpose of this paper is to study the acoustic band structure and transmission spectrum in a periodically 
modulated quasi-one-dimcnsional waveguide, as depicted in Fig. 1. The system has a finite (infinite) extension along 
the y (z) direction and is periodically modulated, along the x direction, by the addition of double stubs, in general 
asymmetric, with different elastic properties than those of the main waveguide. The motivation stems from recent 
studies with interesting results pertinent to electronic JsJ and photonic [^) waveguides modulated in the same fashion. 
Using the transfer-matrix technique we demonstrate the tunability of the acoustic band gaps as a function of various 
parameters of the system, e.g, the length and/or width of the stubs. 

The rest of the paper is organized as follows. In Sec. II we introduce the formalism for studying the wave propagation 
in the allowed polarization and present the necessary details of the transfer-matrix method ]To| . In Sec. Ill we present 
several illustrative numerical results on the band structure and transmission spectrum under various material and 
geometric conditions. The final section contains the concluding remarks. 

II. FORMALISM 

This section is divided into two parts. First we embark on the polarization pertinent to wave propagation in a 
two-dimensional (2D) system. Then we present the transfer-matrix method for quite a general geometry of the unit 
cell of acoustic stub tuners shown in Fig. 2. 

A. Polarization of the wave 

We start with the general equation of propagation of harmonic acoustic waves in an isotropic three-dimensional 
(3D) homogeneous medium 

(\ + p)V(V-u) + pV 2 u + pu} 2 u = 0, (1) 

where p is the mass density, u> the angular frequency, and A and p the Lame coefficients. The longitudinal, vi, and 
transverse, v t , speed of sound are defined in terms of the Lame coefficients: vi — \/{X + 2p)/p and v t = \J p/p. For 
a 2D system at hand the displacement vector u is independent of the z coordinate and one can take d z = d/dz = 0. 
Then Eq. (1) can be written as 

(A + p)V p (d x u x + dyUy) + pV 2 (u x i + u y j + u z k) + pu?(u x i + u y j + u z k) = 0, (2) 

where V p = id x + jd yi and i, j, and k are the unit vectors along the x, y, and z axes respectively. This equation is 
equivalent to the independent equations 

(V 2 + k 2 )u z = 0, (3) 

with k t — Lu/v tl and 

(A + /i)Vp(Vp-Wp) + p^pU p + puj 2 Up = 0. (4) 
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Here the subscript p is assigned to the quantities which qualify only in the x-y plane. It is thus quite reasonably 
understandable that a 2D system can support two independent modes: the out-of-plane modes and the in-plane 
modes, described, respectively, by Eqs. (3) and (4). Equation (4) can be further simplified as follows. We write 

u p = V p <j) + V p x i/j (5) 

with ijj = (0,0, ip). Then Eq. (4) further splits into the equations 

(V£ + kf)</> = 0, (6) 

with ki = lo/vi, and 

(V* + fc t 2 )V = 0. (7) 

Interestingly, Eq. (7), which describes the transverse in-plane vibrations, has formally the same structure as the one 
for the out-of-plane vibrations, Eq. (3). Also, Eqs. (3), (6) and (7) are formally identical to the scalar wave equation 
for the TE polarization in photonic crystals ||. 

It should be pointed out that the splitting of Eq. (4) into Eqs. (6) and (7), is valid only for a homogeneous medium. 
For an inhomogeneous medium, when A and /i are functions of position, this is no longer possible If the system 
is piecewise homogeneous, as the one we consider in Sec. Ill B, problems arise in the application of the boundary 
conditions at the interface of different regions fl2|| that make the separation of the in-plane modes in pure longitudinal 
and transverse ones impossible. 



B. Transfer-matrix technique 



For the sake of generality, we start with a crossbar- like geometry of a single unit cell, as shown in Fig. 2. The origin 
of the Cartesian coordinates is at the uniaxial line intersecting perpendicularly the left arm of the stub of width b and 
length h. The center of the asymmetric stub lies at (x — 6/2, d). We denote the width of the left (right) waveguide 
segments by c (a) and take the x axis parallel to the direction of propagation. We are interested in the solution of 
the wave equation for the out-of-plane vibrations in the form 

V 2 ( j) + k 2 (j)^0, (8) 

where <f> = u z , k = kt, and V 2 = V 2 . It is very important to note that we consider, for the sake of simplicity, 
that the outer medium containing the said acoustic device is made up of some high-density, infinitely rigid material. 
The resulting situation is equivalent to that attained in the case of similar electronic devices surrounded by infinitely 
repulsive walls || . 

In order to solve the scalar equation and describe the system we use the same transfer-matrix method that was 
employed in the study of electronic [fTo[|s|| and photonic |)| tuners. The method relates the incoming to the outgoing 
wave across the stub for arbitrary initial conditions. Inside the waveguide segments, since the solution must vanish 
on the walls, the y dependence is sm(mr(y + c/2)/c) for the left segment, for example. Here the integer n defines the 
number of modes in the respective waveguide. When the two segments connected with the stub have different widths 
and elastic properties, the respective solutions are given by 

<h = V^ 0mX + ci»e" l/3m 1 sin (c m (y + c/2)) (9) 

on the left segment and by 

sin(a„(y + a/2)) (10) 
on the right segment. Here a n — nix j a. c m — mn/c, 



a n = \jk\- a 2 n , ki=w/vi (11) 

and 



/ \ ^ „ ia n (x~b) i — — ia„(x — b) 
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An = \Jkl-c 2 m , k 2 = w/v 2 , (12) 

with v\ (V2) the transverse speed of sound of the material in the right (left ) segment. Inside the stub, </> must vanish 
at y = d — h/2 and y = d + h/2; thus the basic y dependence is sin (kiry±/h), where y± = y T h/2 — d. However, 
the internal solution should also vanish at each side of the stub outside the main segments, and smoothly connect to 
the external one across the contact boundaries between the stub and the segments (at x = and x = b). We first 
construct two auxiliary sets of solutions to the wave equation, one of which matches the waveguide on the left and 
one on the right, with each vanishing elsewhere on the boundary. The appropriate boundary conditions are 

XkL (x = b,y)=0 (13) 

0, y > c/2 

XkL (x = Q,y)={ sm(c k (y + c/2)), -c/2 < y < c/2 (14) 

0, y < -c/2 



and 



x kR (x = o,y) = o (is) 

0, y > a/2 

X kR (x = b,y)={ sm(a k (y + a/2)), -a/2 < y < a/2 (16) 

0, y < -a/2, 



where the subscript k (j = L,R) refers to the number of modes (left/right segments). The solutions x kB are expanded 
as 

X kR = ^2 K sin(7„x) + v n cos(7„x)] sin (/i„y_) , (17) 

n 

where y_ = y + h/2 — d, h n = nn/h, and 



In = V /k s - h n . k s =L)/v s , (18) 

with v s the transverse speed of sound of the material in the stub. The boundary condition at x = requires v n = 0, 
whereas the condition at x = b yields 

r 0, y > a/2 

^ u n sin(7„6) sin(/i„y_) = < sin(a fe (y + a/2)), -a/2 < y < a/2 (19) 
n I 0, y < -a/2. 



This is a Fourier expansion with the coefficients u n given by 

2 

h sin(7 m 6) 

where 

f+a/2 



u m= , L j L . m=l,2,..., (20) 



dy sin(a k (y + a/2)) sin (a m2 /_) (21) 

-a/2 



/2 

1 

[sin (kir - h m s + ) + sin (h m s^)] 



2ir \(k — ma/h) 
1 



(k + ma/h) 

and s± = h/2 ± a/2 — d; then from Eq. (15) we have 



[sin (kir + h m s + ) - sin (h m S-)} > (22) 



2 ^ sin(7 m a;) 

Following the same procedure, we find that 



^-^E ^-5) 7 - sin ^-) ( 23 ) 
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2 ^ sin(7 m (6 - x)) L . 
XkL ^h^ sin( 7m 6) 7 ^ Sln ^-) 



(24) 



where lfc m is defined just as with a replaced by c. The actual wave function in the stub region can be expanded 
in terms of these auxiliary solutions x k R an d Xhl 



— ^y^j {fkXkL + fkXkn) ■ 



(25) 



The continuity of the wave function at x = and b requires fk = C\ k + C\k and f k = c 2 k + c 2 k- Thus one can write 
2 



h 



E 

km 



(C 2 fe + C 2 fc) JT I km + (Clfe + Clfe) — - h 



sin(7 m 6) 



sin(7 m &) 



km 



sin (h m y-) . 



Similarly, matching the derivative at x = gives 
E( Cln ~ sin ( c "( y + c / 2 )) 



2 (c 2fe + c 2k )Ik m - (cu + ci fc ) cos(7 m &)/ ( 



km 



km 



sin(7 m 6) 



7 m sin (/i m y_) , 



and at x = b, 

E( C2 « - c 2 „) ia n sin (a„(y + a/2)) 



2 (c 2 fe + c 2fe ) cos(7 m 6)/ fc R m - (cife + cife)7; 



h ^ 



km 



km 



sin(7 m 6) 



7 m sin(/i m y_) . 



Multiplying Eq. (27) by sin(c/(y + c/2)) on both sides and integrating from — c/2 to c/2 gives 

ch 



, - x. fl 4 V- ( c 2fc + cafc)/)^ - (cut + cu) cos( 7m 6) j£ m L 
1 1Z 1/j ^ ~ rh 2^ sin( 7m 6) 7m /m ' 



km 



Similarly, multiplying Eq. (28) by sin(a/(y + a/2)) on both sides and integrating from —a/2 to a/2 yields 



( C2i -c 2 ^=-x: 



4 v-^ (c 2fe + c 2 fe) cos(7 m 6)/^„ - (cu + cife)J ( 



/cm ri? 



km 



sin(7 m &) 



ImJ-lm- 



We define 



Cfe = Cfc ± Cfc 



to cast Eqs. (29) and (30) in the form 



W c u = 7h E [4A ~ c+ cos( 7m &)4 L m ] 



and 



c/l z_, sin(7m6) 



We now define matrices A, B, D, E, and a whose elements are 



all _ 4 V- cos(7mfr) , IfI 
< fe ~ cfe ^ sin( 7m 6) 7m ^»^» ' 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 
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a ik - — 



ch sin(7 m 6) 



T L jR 



(35) 



nRL — ^ \ " tR t 



L 

km "> 



(36) 



t?RR 



cos(7 m 6) 



ah^-f sin(7 m 6) 



<m 1 lm 1 kr. 



Im km ' 



(37) 



aik = iaiSik , Pik = 0iSik- (38) 
We also define the column vectors C 4 + and whose elements are cf k and c~ k , respectively. In this notation, we have 

0Ci = -AC+ + BC+, (39) 



aC 2 " = -DC+ +EC+, 



(40) 



where A, B, D, and E are real- valued matrices. These two equations determine and C 1 in terms of C 2 and C 2 . 
The result is 



where M^" is the resulting transfer matrix with matrix elements 



M 



/3a 
12 



(41) 



(42) 



M. 



pa 
21 



-$- 1 ad~ 1 e + 



M 



0a 
22 



(43) 



This transfer matrix relates the incoming to the outgoing wave across the stub for arbitrary initial conditions. 

The other building block of a multiple-stub system is a stublcss waveguide segment. The transfer matrix induced 
by the segment of length lij connecting the ith and jth stub is a special case of the matrix M^ a in Eq. (41). This is 
obtained by considering the special case c = a = h, d = 0, a n = /?„, and k 2 = k\ = k s . The result is 



P(ij) 



cos(alij) —ism(alij) 
-i sm(alij ) cos(a/jj ) 



where j = i ± 1. Given M$ a and P, the total transfer matrix for an n-stub system is 



n-l 

M T = Yl \M 0a {i)P{i,i + 1)] M 0a (n) 

i=l 



(44) 



(45) 



Suppose the incident wave, at the entrance of the first stub, is represented by {C^ n , C in } and the outgoing one, at the 
exit of the last stub, by {C^ ut , C~ ut }. Then 



ci 



M 1 



c: 



(46) 



We now discuss the physical conditions imposed on the incoming and outgoing wave components. Depending on 
the frequency, a n and /?,, could be real or pure imaginary. We choose them to be positive for an open channel and 
lying on the positive imaginary axis for a closed channel, so a complete set of eigenfunctions of Eq. (8) is obtained. 
With this choice, C „ ti „ represents either a leftward moving wave or an exponentially divergent wave at positive 
infinity. Now, for the sake of the argument, we will take the incident wave to come from the left; in this situation, 



we must set C„ 



and thus C^ ut 



a 



out 



C out . Physically c out , n represents the transmitted wave components in 
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the nth mode. In the entrance region the components Ci n>n are allowed. Here Ci n>n represents the amplitude of the 
reflected wave in the nth mode in the open channel. For a closed channel the reflected wave is a transient and decays 
exponentially. Explicitly, 

G^n = Gin + Cin = ^l\\G ut + M^Cout, (47) 
Gin = Gin — Gin = M^Cout + M^Cout, (48) 

where are blockwise submatrices of M T . Adding these two equations we can determine the transmitted amplitudes 
G ou t by solving 

2Cin= fe^J Cout ' (49) 
Following this, the reflection coefficients are given by 



l - (MS + M£ - Ml - Af 2 T 2 ) C out . (50) 



The total transmission and reflection coefficients are then given by 



and 



i?= ^" £{ — >' m '" 5 r Q " . (52) 

For the sake of completeness, it is noteworthy that for the bound-state calculation, one must solve the homogeneous 
version of Eq. (49). 

We also stress that in this formalism, once we know how to handle the single stub case, the multiple-stub problem 
can be dealt with little additional labor. This is not the case with either the recursive Green-function method or the 
usual mode-matching approach. 

Finally, the band structure for a periodic system of acoustic stub tuners is computed by solving the standard 
eigenvalue equation 

M B <& = e ik * L I§, (53) 

where k x is the Bloch vector and L = b + 1 is the period of the ID superlattice of acoustic stub tuners. The matrix 
M B is the same as M T for the unit cell but with c = a and f3 n =a n , / is the unit matrix of the same order as M B , 
and $ is the column eigenvector. So the strategy of the computation is that we input the dimensionless frequency 
Q, = u>L/irv wg and calculate w = exp (ik x L), yielding k x — —(i/L) ln(ui). For |iy| = 1 1) one obtains bands (gaps) 
in the band structure for a given set of material and geometrical parameters. It should be pointed out that in order 
to make this paper as self-consistent as possible we have heavily relied on Ref. mn in this section. 



III. ILLUSTRATIVE EXAMPLES 



For the sake of clarity we discuss the numerical results in two parts. First, we consider the case when the waveguide 
and the stubs are made up of the same material. Clearly the band structure and/or transmission spectrum in this 
case reveals the influence of the various parameters involved in the problem. Then we take up the case when the 
materials in the waveguide and the stubs are different. Practically speaking, this case is more complex but richer 
than the previous one in the sense that one has more options to modulate the band structure and/or transmission 
spectrum. We have chosen carbon and epoxy resin as the suitable materials the acoustic system considered is made 
of. This is because these are the materials whose combination was first demonstrated to give rise to a complete band 
gaps, i.e., independent of the polarization of the wave and of the direction of propagation, in 2D periodic phononic 
crystals ]l3]| . The parameters used are p — 1.75 (1.2) g/cm 3 and v t — 711,095 (115,830) cm/sec for carbon (epoxy). 
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A. Same material in waveguide and stubs 



The left part of Fig. 3 shows the first nine bands for a symmetrically stubbed system made up of epoxy, with 
parameters = a/L = 0.5114, = b/L = 0.4886, hi = h/L = 1.125, and = d/L = 0.0. As one can see, all 
the nine bands are separated from each other by stop bands, or gaps, within which the acoustic wave propagation 
is forbidden. Unlike the other 2D and 3D periodic systems ^|JJ], there is a complete gap below a cutoff frequency 
51 c ~ 1.6 down to O = 0. It is found that this (the lowest) gap persists independent of the values of the variable 
parameters. The existence of all nine gaps is well corroborated by the energy dependence of the transmission coefficient 
for n s tub — 50 on the right part of Fig. 3. The numerical results clearly reveal the zeros and ones in the transmission. 
It is noteworthy that the band structure in this figure contains both direct and indirect gaps. For instance, the second, 
third, fifth, sixth, and ninth gaps are direct, while the rest are indirect. We consider it more appropriate to include 
the lowest (and also the widest) gap in the category of direct gaps. The most important aspect of these results is the 
cutoff frequency f2 c below which no propagation at all is allowed. 

Figure 4 depicts the dimensionless gap widths of the three lowest gaps of Fig. 3 as a function of the dimensionless 
stub width bh and the stubs are symmetric, i.e., d — 0. The width of the lowest gap Ai decreases gradually, by 
approximately 45%, with increasing b^ but still remains finite for b — > L. We notice that the maximum of Ai is 
near the cutoff frequency of the waveguide segments A w — 1/cll = 1-96 and the minimum of Ai is near the cutoff 
frequency of the stub segments A s = l/h^ = 0.89. In general, the cutoff frequency of the combined system is between 
the one of the waveguide segments and that for the stub segments. The second lowest gap A2 reaches a maximum for 
bh — 0.7 and decreases slightly before approaching the final minimum at b — ► L. As can be seen, this gap increases 
enormously relative to its value for zero stub width. Similarly, the third lowest gap A3 reaches one maximum at 
bi ~ 0.44, it then vanishes at &l — 0.68, and finally reaches a maximum at b^ — > 1. Note that all three gaps start 
opening up at a vanishingly small but finite value of &l and that the lowest gap remains the widest one over the whole 
range of the stub width. 

Figure 5 represents the band structure and the corresponding transmission spectrum for asymmetric stubs, with 
d = 0.25. The rest of the parameters are the same as in Fig. 3. We observe that the asymmetry has introduced two 
important effects. First, the number of bands accommodated within the same frequency range has increased, from 
nine to fourteen. Second, the band width of most of the bands has reduced. Overall, the effect of the asymmetry 
seems to result in gaps larger in number but shorter in width. This is true despite some exceptions, for instance, the 
case of the lowest gap, which now extends from fl — to ~ 1.7, instead of up to f2 ~ 1.61 in Fig. 3. All gaps in the 
band structure ( left part of Fig. 5) are seen to be well substantiated by those in the transmission spectrum (right 
part of Fig. 5) for n = 50. 

As a function of the asymmetry parameter d the three lowest gaps vary very little, by at most 10%, for < d < 0.3. 
Their dependence on bi is similar to that shown in Fig. 4. As a function of the stub length Hl their behavior, shown 
in Fig. 6, is similar to that in Fig. 4 for Ai and A 2 but somewhat different for A 3 . 

Figure 7 shows the transmission spectrum versus reduced frequency for a system made up of eleven symmetric stubs 
with the central (sixth) stub longer ( = 1.395) and wider (&l = 0.6305) than the rest of the identical, in width and 
height, stubs. When not identical to the other stubs, this central stub constitutes a defect. The solid (dotted) curves 
correspond to the presence (absence) of this defect. One can see that there are five complete gaps in the spectrum 
within the given frequency range. In addition, there is also a pseudogap, centered at £1 ~ 4.35, that corresponds to 
the low transmission or density of states. The defect introduces sharp transmission peaks, marked by arrows, within 
the first four gaps in plane analogy with the electronic || and photonic || case or with that of surface modes of a 
truncated superlattice 14 1 . Another interesting consequence of introducing a defect in the system is the appearance 
of antiresonances |J, such as the one appearing in the fourth band at SI ~ 3.92. We have noted similar effects in the 
case of an asymmetric (d 7^ 0) defect introduced in the system. 



B. Different materials in waveguide and stubs 



We now present numerical results for a system in which the waveguide is made of epoxy and the stubs of carbon. 
Figure 8 shows the band structure and transmission spectrum for symmetric stubs; the parameters are — 0.9, 
bi = 0.9, fiL — 1.5, and d^ — 0.0. We note that there are only seven bands accommodated in the frequency range 
< VI < 20, and every pair of bands has a full gap in between. Moreover, the lowest acoustic gap extends from zero 
to the cutoff frequency fl c ~ 2.3. Some of the bands, such as the fourth and sixth, are seen to be almost flat and 
hence have vanishingly small group velocity. All gaps in the band structure (left panel) correspond well to those in 
the transmission spectrum for n = 50 on the right panel. A prompt comparison of Figs. 8 and 3 reveals that one can 
achieve wider gaps in the band structure if the segments and stubs are made up of different materials. 



8 



Figure 9 shows the three lowest gaps as a function of the stub width bh for the system specified in Fig. 8. There 
are several noteworthy points. First, the lowest gap is the widest one and the third gap is wider than the second 
one until 6l ~ 0.49, where Ai = A 3 . At 6l ~ 0.64, the width of the lowest and second lowest gaps are equal, i.e., 
Ai = A 2 . Also A 2 = A 3 at ~ 0.78 and 0.94. The second lowest gap is the widest one in the range specified by 
0.78 < 6l < 0.94. The third lowest gap vanishes at b^ ~ 0.87 but reappears for 6^ > 0.87 and becomes the widest in 
the range 0.94 < b L < 0.98. Also, Ai = A 3 at b t ~ 0.83 and 0.91. Moreover, Ai = A 2 at b L ~ 0.96 and Ai = A 3 
at 6l ~ 0.98. Finally, Ai reaches a maximum whereas both A 2 and A 3 vanish for b — > L. Similar to Fig. 4, the 
minimum of Ai here is near the cutoff frequency of the waveguide segments A w = 1/cil = 1.1 and the maximum 
of Ai is near the cutoff frequency of the stub segments A s — v a /(v\hL) = 4.1. Again the cutoff frequency of the 
combined system is between the one of the waveguide segments and that for the stub segments. The relative position 
though depends on b^ and the velocity contrast and can be outside the range (A W ,A S ). 

Figure 10 depicts the band structure and transmission spectrum considered in Fig. 8 but with asymmetric stubs. 
The asymmetry parameter is di, — 0.25 and the other parameters the same as those used in Fig. 8. Now there are eight 
bands and eight gaps in the band structure. The lowest gap now extends from O = 0.0 to ft ~ 2.4. The asymmetry is 
seen to have brought about a number of interesting effects. A larger number of bands is accommodated in the same 
frequency range but the bandwidth is reduced. Again, the gaps in the band structure (left panel) correspond well to 
those in the transmission spectrum for a system with n — 50 stubs (right panel). 

As a function of the stub width &l the three lowest gaps behave qualitatively as those of Fig. 9. More important 
is their dependence, shown in Fig. 11 (a), on the velocity contrast between two materials. As can be seen a wide 
modulation can be achieved in this asymmetric structure by just changing the ratio v g /v wgi more than two orders of 
magnitude for A 2 . A less pronounced variation of the same gaps is shown in Fig. 11 (b) as a function of the stub 
length h L . As can be seen, Ai remains almost insensitive to changes in h L and A 2 changes by at most 30%; however, 
A 3 can change by a factor of 10 reaching a maximum at Kl ~ 2.4. 

Figure 12 shows the transmission spectrum versus reduced frequency for a symmetric defect introduced in an 
otherwise periodic system with seven stubs (n = 7). The central (fourth) stub is defect in the sense that its length 
(hL = 3.0) and width (bL — 0.4) are different than those of the rest of the stubs. The other parameters are ol = 0.6, 
bL = 0.2, hL = 1.4, d,L = 0.0. The solid and dotted curves correspond, respectively, to the presence (absence) 
of the defect. There are five complete gaps in this frequency range in the spectrum before introducing the defect. 
Inserting this single defect in the system gives rise to one peak in the third gap and another in the fourth gap. These 
transmission peaks correspond to defect modes similar to those appearing in Fig. 7. 

Finally, Fig. 13 illustrates the transmission spectrum versus reduced frequency f2 in a symmetric system made up of 
one (top panel), two (middle panel), and five (bottom panel) stubs. We remind that the waveguide segments (stubs) 
are made of epoxy (carbon) materials. The parameters are = 6^ = Kl = 0.9. That is, the waveguide has no stubs. 
The lowest panel shows the corresponding band structure. For n = 1 (top panel), the transmission coefficient becomes 
very small but it never approaches zero. In this sense we have only pseudogaps, not full gaps, in the system. As n 
increases the pseudogaps gradually turn into complete gaps (with transmission equal to zero) centered at almost the 
same midgap frequency. It has also been observed that the number of such complete gaps increases with increasing 
n. 



IV. CONCLUDING REMARKS 

We have investigated the existence of tunability of complete spectral gaps in the band structure of a quasi-one- 
dimensional waveguide with double stubs periodically grafted at N equidistant sites. The waveguide segments and 
the stubs can be made up of same or different materials. The latter case is found to be relatively more interesting 
since there one has more and better options to achieve the complete gaps. The single symmetric defect is shown 
to introduce extra modes in the gaps of an otherwise periodic system that result in transmission resonance peaks. 
While the computation of the band structure requires an infinitely long periodic system, the transmission spectrum 
is calculated only for a finite system. As shown though, the gaps in the band structure correspond well to those in 
the transmission spectrum. Accordingly, we conclude that the transmission spectrum in all cases remains consistent 
for n > 10. 

The numerical results we presented pertain only to the out-of-plane modes. However, as mentioned already in Sec. 
II, A, the equations describing the in-plane modes have the same structure in a homogeneous medium. Explicitly, the 
transverse in-plane modes have exactly the same band structure as the out-of-plane modes. As for the longitudinal 
in-plane modes, their band structure and transmission can be obtained directly from Figs. 3 - 7 by changing only 
the frequency scale when the waveguide and the stubs are made of the same material. If this material is different, the 
separation of the in-plane modes in longitudinal and transverse is no longer possible for the structures of Sec. Ill B 
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(see Sec. II A). 

The version of the transfer-matrix method employed [fiof to accomplish the present investigation is simple and 
efficient [p|-^0[ for solving scalar equations. Its efficiency has already been demonstrated in serially connected electronic 
[pLfloJ, optical ||, or the present acoustic devices. As stressed already in its presentation, the method offers important 
advantages in comparison to other mode-matching or recursive Green's function techniques applied to similar problems. 

In principle, the systems we studied here have potential applications in the designing of transducers and ultrasonic 
filters. We hope that the present findings will be tested in future experiments. 
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FIG. 1. Schematics of a quasi-one-dimensional periodic waveguide. The double stubs can be made of the same or different 
material than that of the main waveguide. L = b + I is the period of the system. 

FIG. 2. Schematics of a general unit cell with asymmetric stubs. 

FIG. 3. Band structure (left panel) and transmission spectrum (right panel) for a system with the same material, epoxy, in 
the waveguide and the stubs. The reduced wave vector and frequency are defined by k x L/n and = LQ/nvi, where vi is the 
transverse speed of sound in the waveguide. Notice the lowest acoustic gap below the cutoff frequency fl c ~ 1.6. For the plot 
of the transmission a system of fifty (n = 50) stubs was considered. 

FIG. 4. The widths of the three lowest gaps as a function of the stub width 6_l = b/L. The solid, dashed, and dotted curves 
refer to the lowest (Ai), second lowest (A 2 ), and third lowest (A3) gaps. The material and the rest of the parameters are the 
same as those in Fig. 3. 

FIG. 5. Same as in Fig. 3, but for asymmetric stubs with asymmetry parameter (Ll = 0.25. Notice the lowest acoustic gap 
below the cutoff frequency f2 c ~ 1.7. 

FIG. 6. The widths of the three lowest gaps as a function of the stub length Kl = h/L. The curves are marked as in Fig. 4. 

FIG. 7. Transmission spectrum for a system with eleven symmetric double stubs (n = 11). The central (sixth) symmetric 
stub is a defect with Kl = 1.395 and &l = 0.6305. The rest of the parameters are the same as in the previous figures. Notice 
that the defect creates new modes or states in the gaps in an otherwise defect-free system. The peaks of these modes are 
marked with arrows. 

FIG. 8. Band structure (left panel) and transmission spectrum (right panel) for system with symmetric stubs with the 
waveguide segments (stubs) made of epoxy (carbon). The parameters used are ol = 0.9, &l = 0.9, and Kl = 1.5. Notice the 
lowest acoustic gap below the cutoff frequency Q c ~ 2.3 and extending down to Q — 0.0. 

FIG. 9. The widths of the three lowest gaps versus the stub width 6l for the system studied in Fig. 8. The solid, dashed, 
and dotted lines refer, respectively, to the lowest (Ai), second lowest (A2), and third lowest (A3) gaps. Notice especially the 
strong variation of A2 and A3. 

FIG. 10. Same as in Fig. 8, but for a system with asymmetric stubs ( (Ll = 0.25). The lowest acoustic gap occurs below the 
cutoff frequency f2 c ~ 2.4 and down to = 0.0. 

FIG. 11. (a) The widths of the three lowest gaps as a function of the velocities' ratio v s /v wg . (b) The three lowest gaps as 
a function of the stub length Kl- Notice the strong variation of A2 and A3, especially in (a). 

FIG. 12. Transmission spectrum for a system with seven symmetric (d — 0) double stubs (n = 7). The waveguide segments 
(stubs) are made of epoxy (carbon). The central (fourth ) stub is a defect with Kl = 3.0 and 6_l = 0.4. The other parameters 
are a,L = 0.6, 6l = 0.2, and Kl = 1.4. The solid (dotted) lines refer to the transmission with (without) the defect. The defect 
gives rise to modes in the gaps in an otherwise defect-free system. The peaks of these modes are marked with arrows. 

FIG. 13. Evolution of the transmission spectrum as a function of the number of stubs n for a system with waveguide 
segments (stubs) made of epoxy (carbon). The parameters are ol = £>l = tiL = 0.9. Notice that as n increases the pseudogaps 
gradually turn into sharply defined complete gaps. The lowest panel shows the band structure. Notice that the first band is 
very narrow. 
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